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In the context of a family of scalar-tensor theories with a dynamical A, that is a binomial on 
the scalar field, the cosmological equations are considered. A general barotropic state equation 
p — (7 — l)p, for a perfect fluid is used for the matter content of the Universe. Some Friedmann- 
Robertson- Walker exact solutions are found, they have scale factor which shows exponential or power 
law dependence on time. For some models the singularity can be avoided. Cosmological parameters 
as Sim, f^A, qo and to are obtained and compared with observational data. 



I. INTRODUCTION 

Unification theories have a nonzero cosmological constant that is about 120 orders of magnitude larger than the 
observed value for A, this constitutes the cosmological constant problem [Q- In order to explain and solve such 
a problem, and to make compatible the actual observational data with the inflationary scenario and particle physics 
expectations, a time dependent cosmological constant was proposed This old idea has received a lot of attention 
(see e.g. [Q- |4^ ). What people have in mind is to make the vacuum energy dynamical. In such a way, during the 
evolution of the universe, the energy density of the vacuum decays into particles, thus leading to the decrease of the 
cosmological constant, obtaining as a result, although small, a creation of particles. 

A broad summary of cosmological models with time dependent cosmological "constant" is given by Overduin and 
Cooperstock [Q, re-examining there, the evolution of the scale factor when A is given as function of t, a{t), H, or 
q. A fairly general equation of state is considered and new numerical solutions are obtained, but as in most of the 
previous works, the time dependence of the cosmological term is introduced ad hoc. 

An alternative is an effective time dependent cosmological "constant" in the context of a scalar-tensor theories, 
which becomes a true constant for i ^ [Q. Using Jordan-Brans-Dicke theory (JBD) in particular, the "graceful 
exit" problem of old inflationary cosmology might be improved. It remains the problem of determining the JBD 
parameter w, that according to solar system experiments its value is ||u;|| « 500, which has been derived from timing 
experiments using the Viking space probe . A better estimation of this parameter should be obtained from measure 
of others cosmological parameters in order to constrain to more strongly than by means of solar system experiments 
^7\ . However, theories of very early universe as string theory, are better described in the context of JBD, and shows 
that uj can take negative values |Q. 

Thus, scalar-tensor theories, and in particular JBD, are better theories, in order to get, in a natural way, a time 
dependent cosmological constant. Clearly, recent observational results restrict this kind of theory, e.g. the type la 
super nova (SN la) results, which in 1998 show that fl\ ^ 0.6 implying that our universe is speeding up. Thus, 
a model which attempts to describe the cosmological constant behavior, should take into account the observational 
evidences. 

In a recent work we investigated the effect of a time dependent cosmological constant, in a family of scalar-tensor 
theories. There, we get cosmological models in the coasting period, where the time dependence on the cosmological 
constant occurs in a natural way. In such a models we assumed a simple relation X{4>) = c(f>{t)"' , (with c and n 
constants) . 

The existence of inflationary phase in scalar-tensor theories (STT) has been investigated by Pimentel and Stein- 
Schabes |5^, finding infiationary phases for a polynomial cosmological constant in a general STT, which includes 
Brans-Dicke model with non-zero cosmological constant. On the other hand, Guendelman has investigated the 
requirements of the potentials in order to have scale invariance. There was found the form of the potential needed 
by the global invariance, which in addition, its energy in the conformal Einstein frame has the characteristics for a 
suitable inflationary universe and A decaying scenario for the late universe. 

Motivated by these ideas, we shall consider a general STT as in our previous work but now we shall consider a 
binomial A function on 0(i), in order to obtain exact solutions of the field equations, from which we obtain some kind 
of inflationary cosmological models and related cosmological parameters. In fact, we obtain in most of our solutions, a 
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power law growth for the cosmofogical scale factor a{t) t" ^ where cr ^ 1 implies inflationary models. As it is known, 
this is a generic feature of a class of models that attempt dynamically to solve the cosmological constant problem. In 
our models cr is a free parameter (at least in most of our models), in order to be adjusted by physical conditions and 
to be in agreement with recent data for SN la, that implies cr « 1, and which is consistent with the nucleosynthesis 

il- 

Most of our solutions predict an accelerated expansion, such solutions are in agreement with the SN la results, but 
rim and f2A, depend on free parameters of our model. In some specific cases we get solutions with exponential growth 
of the scale factor. 

In section || we obtain the field equations and introduce our main ansatz. Section [II is devoted to obtain an 
expression of the density parameter and the corresponding contributions of the matter and scalar field. In section |^ 
we consider the vacuum case and obtain a set of exact solutions as well as cosmological parameters. In section ^ we 
consider the general case of a barotropic equation of state and obtain exact solutions for this general case. As example, 
we calculate specifically the case of a dust fluid and a stiff matter fluid. Also, we get and discuss the solutions of the 
radiation case in section VI and false vacuum case in section VII. Finally in section VIII we summarize our results. 



II. FIELD EQUATIONS 



We start with the action for the most general scalar-tensor theory of gravitation [Sij 
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d^xV^[^R - <f>-^Ljg'"'d^(f>d,^ + 20A(0)] + Sng, 



(2.1) 



where g — det ( g^^), G is Newton's constant, Sng is the action for the non-gravitational matter. We use the 
signature (—,+,+,+). The arbitrary functions Lo{(t>) and A((/)) distinguish the different scalar-tensor theories of 
gravitation; X{(j)) is a potential function and plays the role of a cosmological constant, w(0) is the coupling function 
of the particular theory. 

The explicit field equations are 
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where G^i, is the Einstein tensor. The last equation can be substituted by 
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where T = T^ is the trace of the stress-energy tensor. In a previous work was demonstrated that the divergence- 
less condition of the stress-energy matter t enso r is satisfied if the field equation (2.3) is satisfied too, although our 
field equations are given by Eqs. ( p^ ) and (2.4). 

In what follows we shall assume Lo{(j)) — constant^ A = A(0). The corresponding field equations with a perfect fluid 
for the matter content in the isotropic and homogeneous line element 
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will be considered. Thus the field equations are 
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(3 + 2c^)-2(a-0^) -^(p-3p) = 0. 



Where we have assumed </> = (/>(i), and the derivatives respect t arc denoted by a dot. 

Assuming a barotropic equation of state p ={'-/— l)p and transforming to the time t defined as 

T = [ <j)^/^dt, 



-0, 



the set of equations (p.q)-( |2.§| ) are rewritten in the following way 
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where the derivatives respect to t are denoted by a prime. 
In what follows we shall consider two important assumptions: 
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where m, a, Ai, A2, ni and n2 are constants. The first assumption is a very well known one (see e.g. [p5| and 
reference therein), and it has been used as a condition for the deceleration parameter to be constant for flat models 
in Brans-Dicke theory. Furthermore, with this condition our field equations simplify notoriously and allows us to 
obtain exact solutions. The second condition is the main assumption of the present work which is motivated by the 
cosmological no-hair theorem for scalar-tensor theories | |5l[ |, in order to study infiationary solutions in a theory of 
gravitation with a naturally dynamical cosmological constant. In this work we always work in the Jordan frame, 
where G is variable, however we could make a conformal transformation to the Einstein frame where G is constant 
and we have General Relativity plus a minimally coupled scalar field, then our potential becomes an exponential one, 
i.e., V1+V2 ~ exp(eni0c) + exp(en2(/'c) (where e is a constant and 4'c is a canonically defined scalar field). This is the 
type of potential according to Guendelman ||5^ , that is necesary to have scale invariance in a theory of gravitation 
free of the cosmological constant problem, that is, one with an early expanding phase and a A-decaying for late times. 
Detai ls of th e conformal transformation for STT can be seen in Ref. With these assumptions, from equations 

(p0|)-(pl2|) we get 
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(2.17) 



In the following sections we shall find exact solutions on different cases, as well as cosmological parameters which 
allow us to compare with actual observations of today value of Hubble parameter Hq, the actual value of deceleration 
parameter qo, the density parameter flm as well as the value of the vacuum energy density parameter fl/^. 
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III. THE DENSITY PARAMETER 



Before we compute exact solutions of the set of field equations ( 2.15 ) 
and ri^, accordin g to our proposed model. 
Assuming fc = 0, Eq.(2_^) is written as 



( 2.17 ), we shall get a general equation for 
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Taking into account the proposed relation ( 2.13 ), we get 



LU 1 

Qm? m 



(3.4) 



According to the SN la observations 1 49 1 , the favored value of il^ 
constant 
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0.4 ± 0.1 is given as a constrain to a cosmological 
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this implies JIa ~ 0.85 ± 0.2. It means that the SN la results are sensitive to the acceleration of the expansion, and 
constrain 4f2m/3 — i^A, which corresponds to the acceleration parameter at the median redshift of this objects, z 0.4. 
Then the combination fip = + is constrained by the microwave background radiation (CBR) anisotropy. So 
that, r^o ~ 1 ± 0.2 obtained from COBE and other measurements (see e.g., |5^), and together with D,m 0.4, define 
a concordance region for J^a ~ 0.6, becoming the best fit for the universe model |5^ . 

In what follows, we shall compute both parameters of density in addition to the exact solutions of our field equations, 
in the different cases which we consider in this work. 



IV. THE VACUUM CASE 



Considering a vacuum case {p — 0), we get from Eqs. (2.15)-(2.17) the corresponding set of equations 
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^(3 + 2w) + (^)^3 + 2cj)(i - 3m) - 2Ai0"i-i - 2A2</)"'~^ + 2Aini<?!>"i-i + 

2A2n2(^"^-i = 0. (4.3) 



Naturally for vacuum, the energy density is due to the contribution of the scalar field 0(t): Vl^ — 1. We found exact 
solutions of this set of equations in the following cases: 
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1. k = 



• m = 1/2 This solution is not relevant for our purpose because becomes null. 

• m = 2/3 
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= Ai </)(0~ 


LO 


= -7/3, 



(4.4) 

where = V^ci, ai — a {2ci)~^^'^ , Ai — cl/2, A2 = 0, and ci is an integration constant. This solution 



was written directly in terms of time t, according to (2.9). Here the expansion factor is decaying with the 



time, in conflict with observations, and lo has a negative values. A discussion of the meaning of negative 
values of uj is given in |Q . 

The Ricci scalar for this case is given by the following expression 

5^=y^-^ (4.5) 

where we can see that there is an initial singularity. The today deceleration parameter has a negative value, 
i.e., the present solution is an accelerated cosmological model 

(70 = -4, Ho = -^t^\ (4.6) 

according to this solution, the actual Hubble parameter has negative values, then we conclude that this 
model has not physical meaning. 

m ^ 1/2, 2/3 

On this case we get two families of solutions: 
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and c' is an integration constant. According to the solution of the present case, A is a monomial 
function on (/)(t), although a time decaying one. Clearly the present solution requires ma < in order 
to be expanding. 

On the other hand, the Ricci scalar is given by 

= 12mcr(4TOcr+ l)t"^ (4.9) 

According to this expression, the corresponding solution for the curvature scalar has an initial singu- 
larity. For this model, the present deceleration and Hubble parameter, are 
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qo = -l--^, Ho = -2ma ito)-\ h = -^2m\\a\\. (4.10) 
Zma Hq 

From these last equations we can see that the present model expands with acceleration if ma < —1/2. 
On the other hand, assuming ffo ~ 65±5 km s~-^Mpc~^ |5^, we get to ^ 15.05 ± 1.96 \\ma\\ Gy, then 
the estimated age from this model is small compared with actual accepted values of to- 

For this particular relation between uj and m, we get the following solution 

if)(t) = c'exp[(?!)2t], 
a(t) — 02 exp[— m(/)2i], 

A(^) = XicPit), (4.11) 

where 02 = c'^^^/i, 02 = ac'^™, A2 = 0, Ai = Ac'{'im? — 3to — uj/'l)v'^ and c' is an integration 
constant. In this case we get an inflationary exponential solution provided that m02 < m < 
or 1/2 < TO < 2/3. This model is non-singular, as we can see from the Ricci scalar 

3fi = 24c'm2(6m2-7m + 2). (4.12) 

The present deceleration and Hubble parameters are given by 

qo = -l, i?o = c'to(12to2 - 14m + 4), (4.13) 

thus, this model expands with constant acceleration from a non-singular state, and with constant 
Hubble parameter. 



2. k^Q 

• TO = 1/2 



0(0 = 01 ^-^ 

a{t) — fli t, 

A(0) = Ai0(t), (4.14) 



where 
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Ai — —5-, A2 — 0. 



According to this solution, a{t) grows linearly with the time at a constant rate. The cosmological "constant" 
A decreases with the time. In order to have a real a(t), we must have fc/(3 -I- 2a;) > 0. This is a coasting 
cosmological solution which has initial singularities as it is shown from the corresponding Ricci scalar 

5R= 12(3 + 2w)i-2, (4 15) 

On the other hand, as we have said, the today deceleration parameter becomes null, and the Hubble 
parameter is given by 

go = 0, H = t^\ (4.16) 
then, to = l/Ho ~ 15.05 Gy, in relative agreement with actual observations. 
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V. EXACT SOLUTIONS FOR THE CASE WITH A BAROTROPIC EQUATION OF STATE 



In what follows we shall consider p ^ 0, so that returning to the set of equations ( 2.15| )-(2.17), we get from Eq. 
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Using this expression in equations (2.16) and (2.17), the set of field equations are reduced to the following two field 
equations 
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From equation (5.2) with 7 7^ (the false vacuum case will be consider in section VII ), we have 
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Using this expression and its derivative in Eq. (5.3) we get the following equation 
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In order to solve this differential equation, we shall consider the value m = 1/2, so that equation ( p. 51) is reduced to 
the following one 



37/(/)'\2 37-4 fc 
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From this differential equation we have the two possible cases: k = and fe 7^ 0. The case 7 = 4/3 will be consider 
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1. For fc = we get the following solution which we write in terms of time t as 

a{t) = ait-"^/^, 

A(0) = Al0(^)-'/^ 



(5.7) 



where 
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and Cl is an integration constant. From this solution the expansion condition is cr < 7 > 2/3. 

In this case, the Ricci scalar and deceleration and Hubble parameters are given by the following expressions 
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Causality requires < 7 < 2, so that this is an accelerated model for 7 < 4/3, just at 7 = 4/3, go = 0. On the 
other hand, to ~ ^/Hq ~ 15.05 Gy for 7 ^ 4/3. Then this is a kind of solution where the cosmic expansion is 
driven by the big-bang impulse. 

The energy density parameter and the contribution of the scalar field (/)(t) are given as follow 
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In order to have a positive values of Um, w < 1/2 is required, including negative values of co. On the other hand, 
flm ^ 0.4 ± 0.1, and < 7 < 2, from our equations for ri,„ and fl^, we get a restriction for ui: —0.1 < uj < 1/2. 

2. Considering now the case k ^ from Eq. ([3.6|), we get the following solution in terms of the parameter r: 



for 



> 



where now a 



(I){t) = Cl cosh {(3t) , fc 7^ 0, 7 ^ 0, 4/3 

(37-4) / k 



(5.11) 



4_3-y: P — '''^^2a'^' \l 3+2lu ' ^^"^ "^1 integration constant. According to Eq. (2.13) we get 



a(r) = fli cosh {(3t) , fc ^ 0, 7 7^ 0, 4/3 

The Ricci scalar on this case is given by 

SR = ri(r2 + rg) cosh'' [/3r] - rirj cosh''-^ [/3r]. 



(5.12) 



(5.13) 



where ai — ac^ r\ — 34^2^ ' ^2 = 3(2 + 7), and rg = 37 + 2 + ^u. In order to know the singularities 
of this solution, we calculate the nonzero curvature invariant [pO[, which for this case, are given by 
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R3 = j^Rl 
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then, it is enough to calculate 



Ri = si{s2 - ssf cosh^'"-^ [/3t] + si(s3 + Si)^ cosh^'' [/3r] + 

2S1(S2 - S3)(S3 + S4)cOSh2'"-2 [^t] . 
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where 



4a4(3 + 2w)2' 

53 = k, 

54 = fc(3 + 2uj). 

From equation ( 5.13| ), — > cx) provided that r — s- ±oo, and at le ast o ne exponent is positive, i.e., a > 
or (7 > 2. On the other hand, from the curvature invariant, Eq. (5.15), i?i ^ oo requires that r ±oo, 
and (T > 2 (7 > 2/3), a > 0(7 < 4/3) or u > 1 (7 > 0). So that the solutions of this case, are singular for 
< 7 < 4/3, and r ±00. 



• for 
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t) = ci cos (/5r) , fc 7^ 0, 7 ^ 0, 4/3 



(5.16) 



a(T) = ai cos (/3r) ' , A: ^ 0, 7 7^ 0, 4/3 



(5.17) 
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where a and ai are defined as in the paragraph under equation (5.11) and (5.13). In this case, P 

5? = ri(r3 - ra) CDs'" [/3t] + rira cos'""^ [/3r], 
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where now 



ri 



3ci 



3 + 2uj 
r2 = 3(2-7), 
r3 = (4-37) + 2||3 + 2tj| 



According to (5.14), as in the previous case, we need to calculate Ri only, which for this case is given by 

i?i = si{s2 + 53)^ cos^'"-'^ [/3t] + si(s4 - 53)^ cos^'^ [i3t] + 

2si(s2 + S3)(S4 - S3) cos^-^-^ [/3r], (5.19) 

where 

3 c? 

Sl 



4a4|3 + 2u;|2' 



52 = ||fc||(2--7), 

53 = ||fc||, 

54 = /c||3 + 2tj||. 
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From Eqs. ( 5.18| ) and ( ^.ig| ), 3? ^ oo as well as i?i ao, if [/3r] — > ±(2n + l)7r/2 and at least one exponent 
on the respective expressions of 3? and Ri is negative: cr < 2 (7 < 2/3), a < 1 (7 < 0) or cr < (7 > 4/3). 
Then we have the two ranges of 7 for which the solutions of the present case, are singular: < 7 < 2/3 
and 4/3 < 7 < 2, since furthermore causality requires 7 to be in the interval < 7 < 2. 



On the other hand, from equations (5.1) and ( p.4| ) we get respectively, for both possibilities ^^^^ 

p = pia(r)-37, 



' <OOT^>0 
(5.20) 



Xi^) = Ai,/.(t) It-i + A20(r), (5.21) 

where 




As we can see, X{(j)) remains as a binomial function of (j> ii j ^ 4/3. In order to analyze the behavior of the 
obtained solution in terms of the cosmological time t, we shall give some examples. 



A. Dust fluid 

One interesting application of a barotropic equation of state corresponds to a dust fluid (7 = 1), on that case the 
solution reads as follows 

1. fc = 

a{t) = ai t^, 
A(0) = Ai0(t)i/2, 

p = pia{t)-\ (5.22) 

Where = cj~^, ai = ac^^, Ai = A{2uj — 5)0^^, A2 = 0, pi = i^^cj^^a^^, and ci is an integration constant. 
This is an extended inflationary solution, with a time decaying cosmological constant and initial singularity, as 
it is shown by the corresponding Ricci scalar 

3? = 36i"2 (5.23) 

The expansion takes place with a constant acceleration 

12 2 

Z to -Ho 

With Hq ~ 65 ± 5 km s~^Mpc~^, we obtain to ~ 30.1 Gy, which is too big value compared with the globular 
cluster age. 

The density parameters for a dust fluid are given by 

" - 3 3 ' 

n^ = Y + l- (5-25) 

On this case, as in the general case, in order to have a positive values of flm, it is required that uj < 1/2, 
including negative values. On the other hand, according to observational results, ^ 0.4 ± 0.1 and flj^ ~ 0.6, 
then oj is restricted to be a; ~ 1/5. 
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2. kj^O 

• For 



3+2u 



> 



The solutions in terms of the time t is given by 

p = pia(t)"^, 



(5.26) 
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Pi 



kc^ 3 

47ra^ ' 



and ci is an integration constant. Clearly a{t) and <^(f) must be positive, in order to be physically significant; 

— 1/2 

this requirement restricts the range of values which t can take: t < (f)^ , and ci > as we can see from 
the definition of 4>i . In this case the cosmological term, in spite of being a binomial function on decays 
with the time. 

The corresponding Ricci scalar and curvature invariant show that this solution is singular: 
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According to our analysis of the general solution, for a = 4/(4 — 87), with 7 = 1 we g et a ~ 4. T aking into 
account the general equations of the Ricci scalar and curvature invariant, Eqs. ( 5.13 ) and ( ^.15 ), 3? — > 00 
and Ri — > 00 for a — A and r — > ±00, which corresponds to a finite value of t, according to the time 



dependent solution (5.26). Furthermore 7 = 1 < 4/3, is consistent with our singularity requirement in our 



discussion for the general case with 



> 0. 



The today values of the deceleration and Hubble parameters are given by the following expressions 



<lo 



-( — 

2 V01 1 



Ho 



201 ^0 
01 to - 1 ^ 




Because 0iio < 1, then qq > 0. This model expands from t ^ —c, 



-1/2 



(5.29) 



until t — 0, then it contracts until 



, in both cases with positive deceleration parameter. The numerical value for to depends on the 



t = 

values of the free constants. 



• For 



3+2uj 



< 



The corresponding solution in terms of the time t is given as 

0(t) =Cl[l + 0l t^^, 
a{t) =ai[H-0i t% 
A(0) = Ai0(O+A20(t)i/2 
P Pia{ty^, 



(5.30) 



where 



— II — 

4a2ll3- 



2a;' 



fli — ajc^ 



1/2 



_ 2fc 
Ai — —5-, 



\i = 



kc^ 



1/2 



kc 



1/2 



Pi 



Attq? 
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and ci is an integration constant. In this case we have not restrictions on the values which t can take. If 
— 1/01 < {to + c)'^, the expansion takes place with non-constant acceleration and without singularity, as it 
is shown from the corresponding Ricci scalar and curvature invariant 



3? 



"•l + 201 



(5.31) 



1 + 01 t2 



k 1 2 

20? + 2(1)1-— . 



(5.32) 



According to our general analysis of this case, from ( 5.18| ) and ( ^.19 ), 3? and Ri do not diverge for cr = 4 
and [/3t] ±(2n-|-l)7r/2, in agreement with the time dependent solution ( 5.30 ) which has not singularities. 
7 = l<4/3in this model, which is consistent with our condition 7 > 4/3 as requirement for existence of 
singularities. 

For this solution, the corresponding present values of the density and Hubble parameters are given as 



201 to 



± 



'^"^ 2 201 *""' ^""0it2 + i' ^"-Ho^MH? 



(5.33) 







As in the previous case, the values of Iq depends on the free constants of our model, but in this case, for 
Hi ~ 01 =^ <o ~ l/Ho. 

B. Stiff matter fluid 



Another interesti ng a p plicat i on of a barot ropic e quat ion of state is a stiff matter fluid for which 7 = 2. In such a 
case, the solutions ([3.7|), (5.11),( 5TT^ ), ( ^.16 ), and ( ^.17 ) arc the following ones 



1. k=0 

The solution in terms of the physical time t is given by 

0(t) =0i(c_t)-i, 
a{t) = ai(c- t)i/^ 
A(0) = Ai0(t)2, 
p pia(t)"^. 



(5.34) 



where now 



„l/2 

~2~' 



aV2 . 4 1 l-2u 

ci Stt ci 



This solution has a physical meaning for t < c, where c is an integration constant and the scale factor increases 
very slowly with the time, and with constant deceleration. The Ricci scalar and Hubble parameter, are given 
by 



= 0, 



(5.35) 



qo = 1, 



1 1 



2 to - c 



(5.36) 



Naturally this model is not valid today because it shrink for the allowed range of t. 
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The corresponding density parameters are 



" - 3 3 ' 

n^ = l + Y- (5-37) 



In order to have a positive values of flm, then uj < 1/2, as we have claimed in the discussion of the general 
solution. The observational results — 0.4 and fl\ — 0.6, determine uj = —0.1. 



2. k^O 

• For 



3+2u 



> 



The solution of this case, in terms of the time t is given as 



(f>{t) = Cl 



1 + 01 r 



1 



a(t) = fli 

A(0) = Ai0(t), 
p = pia(i)~^. 



1/2 



(5.38) 



where 



3+2uj ' 



6i > for Cl > 0), ai 



.-1/2 



, Ai 



2fe 



A, 



0, Pi 



^ a^, and ci is an 



integration constant. Here a{t) grows with the time, from a minimum radius oi there is a non-singular 
state. The expansion takes place with non-constant acceleration, as we can see from the corresponding 
Ricci scalar and curvature invariant, which are given by 



3? = 




(5.39) 



Ri 



Cifc 



kci 



(5.40) 



In this case 7 = 2 implies a — —2. from Eqs. ( 5.13 ) and ( 5.15 ) with this value of cr, neither 3? nor R\ 
diverge for r — > ±00 which is consistent with our require ment 7 > 4/3 fo r the avoidance of singularities. 
This analysis is in agreement with the inspection of Eqs. (5.39) and (5.4C). 
The deceleration and Hubble parameters are 



90 



1 

= -— t, 



' 



Cl 



a2 3 
1 



2uj 1 



Cl 



+2 

q2 3+2w''0 



2iJ, 



± 







a2 3 -f 2w 



ml 



Cl 



k 



(5.41) 



• For 



3+2w 



< 



m 


= Cl 


1-01 


a{t) 


= ai 


1-01 




= Ai0(<), 


p 


= pia(i)-6, 



1/2 



(5.42) 



13 



where 



In order to have a physically solution, it is required that 0i < 



, q2 II 3+2w I 

corresponding Ricci scalar and curvature invariant in the present case, are given by 



c < 0. The 



6 



1 - 01 i2 



fcci 



(5.43) 



Ri 



cik 



kci 



(5.44) 



In this case 7 = 2, then c = — 2; such that in the general solution ( Eqs. ( p. 16 )- ( 5.19| )), — > 00 and 
Ri — > CX3, i.e., the solution for this case is singular provided that [/3r] => ±(2ri + l)7r/2, which corresponds 

According to our singularity discussion under equation (p.l9|), for 7 — 2 the corresponding 



to t 



solution should being singular, as we can verify from inspection of Eq. ( ^.43|) and (5.44) 
the present deceleration and Hubble parameters are 



90 
Ho 

to 



'-0 ' 



n/2 



k 



3 + 2uj 



tr 



ct2 II 3+2w H o 



2Hn 



± 



a2 3 + 2w 



4H§ 



k 



(5.45) 



As we can see, (/>! < =^> < 0, then this model is accelerated. The numerical values of and 
depend on the values of the free constants. 



VI. THE RADIATION CASE 



We shall consider the radiation case for which 7 = 4/3, so that returning to Eqs. ( 2.15 )-( 2.17 ) and following a 
similar procedure as in section^, we get Eq. (5.5) with 7 = 4/3. In order to solve this differential equation, we shall 
assume m = 1/2, then we have for this case 



= 0, 



(6.1) 



for which uj ^ —3/2. This differential equation has the following solution 



(6.2) 



where c and ci are integration constants. According to Equations ( 2.13 ), ( p.l[ ) and (5^) with 7 — 4/3 and m = 1/2 
we get respectively 



a(T) = aie-=^/2^ 

A(0) = Ai(^(t), A2=0, 



(6.3) 



where oi = ac^ Ai 
is given as follows 



-i3 + 2uj) 



3 k 
2^ 



and pi 



3 

IStt 



k _ ujc^ _ 3cj 
2 4 



a . In terms of the time i, this solution 







ait) 


= ai t, 


A(0) 


= ^im. 


P 


= Pia{t)-^ 



(6.4) 



14 



where (j)i = a2 = ^ and A2 = 0. This is a singular solution according to the Ricci scalar which for this case is 
given by 



+ (6^5) 



The deceleration parameter becomes null, instead the today Hubble parameter is given as 

qo = 0, Ho = ^, (6.6) 

and therefore to = 1/Hq ^ 15.05 Gy. 

The case A: = is not excluded from the solution (3.4). For A: = we get the density parameter from the matter and 
the scalar field as 

o 3 

"m — — UJ — — , 

^4,^ uj+^. (6.7) 
It is required that uj < —3/2, in order to have flm > 0, and the observational accepted values of Qm and fi^, determine 



VII. THE FALSE VACUUM CASE 



We analyze now the case 7 = 0. We shall follow a similar procedure as in section then we get from Eqs. ( |5.2| ) 
and (5.3) with 7 = 0, the following set of equations 



— 2m — UJ — 



2k 



2m -1 



= 0, 



(7.1) 



-(3 + 2.)+(^) 



12to + 3m — Quim + 3ui 



1^0-- +2M + 2^ = 0. 
a 6 d(b 



(7.2) 



We consider first the case /c = 



The set of equations (|7.lD-(|7.2|) with fc = has the two possible set of solution depending on the relation between 
TO and Lo: 



1. UJ = —1 — m 



cj){t) = ciexp[0i t], 
a{t) — aiexp[— TO01 1], 

ci 



(7.3) 



where 



„l/2 



(2to— 1), ai=aci™, Xi — ci, A2 



ci 



(2to- l)^(3m- 1). 



This is an inflationary solution if m(j>i < 0. This condition means < m < 1/2, which implies a condition 
on the range of the values of w: 1/2 < a; < 1. In order to have physical solutions, ci > 0, which means 
p < 0. Of course, these solutions have not singularities, as we can see from the Ricci scalar 
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3?= 12cim2(2m- 1)2. 
The deceleration and Hubble parameters are given by 



90 



-1, Ho = m{l - 2m)cy^ 



(7.4) 



(7.5) 



thus the model is accelerated. On the other hand, the density parameter due to the matter and scalar field, 
are 



1 



1 



C2 4)i 2 exp[-4>it], 



c2(/)i exp[-0ii] + 1. 



(7.6) 



Here, it is required C2/C3 < 0, in order to have a positive values of fi^- As we have seen, < to < 1/2 =; 
TO > and < 0, which means fim and increase exponentially with the time, keeping ri,„ + il^ = 1. 

2. 7^ — 1 — TO 

with this condition and fc = 0, the corresponding exact solutions to the set of equations (7.1)-(7.2) are 

0(t) = ic~tr, 

a{t) = ai (c-i)-""", 
A(0) = Ai 0(t)"i + 



where 



1 - 2to 

TO + ti) + 1 ' 



2lu + 3 

n-i — 

2to - 1 



Ai 
A2 



2 
a 



(2to - 1)2 



2TO + tj + 1/2 



(3 + 2tS)(2m + w + 1/2) + (2to - 1)(-I2m2 + 3to - 6tocj + 3tJ + 3/2) 



C2- 



This solution corresponds to an extended inflationary model if ma < 0. In order to have a physical solution, 
it is required t < c. The cosmological constant is a decaying function of time: A ~ + t~'^ , with a > 0. 
But according to the above mentioned condition ma < 0, it is clear that to should be negative. 
This set of solutions is singular just at i = c, as we can see from the Ricci scalar 



3fi — 6to(2to — 1 



4to — 3to — w — 1 



1 



(to + w + i)2 {c-ty 



The present deceleration and Hubble parameters for this model, are given as 

TO + cij + 1 m(l — 2m) , m(l — 2m) 1 

go = 1 - Hq = -—^ — — r(c-^o) , to ^ c 



m{2m — 1) 



TO + CJ + 1 



TO + CJ + 1 Ho 



(7.9) 



According to this results, the expansion of this model takes place in accelerated way if go = 1 + i/ma < 0, 
but we have seen that ma < 0, then ||TOtT|| < 1 is the required condition for accelerated expansion. The 
value of to clearly depend on the value of to and w, but we know that, because the restriction on t, this 
model is not relevant at present times. 
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The corresponding density parameters due to the matter and scalar field for the present case, are given by 

1 1 ^-2/- 



n, 



3m2 a2'^i 
1 1 



Sirpi - C2<p{t) 1-2™ 
Ai + X2(t>it) 1-2™ 



1 

6rn2 777, 



o 2 2 ''^1 

Both values depend on the free constants, but for all time, it is satisfied that fim + = 1- 



(7.10) 



In order to solve the set of equations (7.1 )- ( |7.2| ) with /c 7^ 0, we shall assume m = 1/2, then we have respectively 

Vn2 4fc 



(7) (3 + 2.)--,=0, 



(7.11) 



(7.12) 



From Eq. ( 7.11 ) we get 



solution to this equation is given by 



FromEq. ( |7l3| ) we get 



4fc 



^(3 + 2cj) 



1/2 



2. / k 

(r) = cie° V 3+2c.^ 



4fc 



a2(3 + 2a;)' 



We use this last equation in Eq. ( 7.12 ) from which we get a differential equation for A(0): 



(7.13) 



(7.14) 



(7.15) 



(7.16) 



the solution of this equation reads as follows 



A(<^) = ^0(r) + -^. 

(p(T) 



From Eqs. (2.13) and (5.1) with m = 1/2, we get respectively 



i(t) = aie ° VVk^^^ 



ci 



(7.17) 



(7.18) 
(7.19) 



where oi = aci ^^"^ . In terms of i, this solution becomes 



17 



a{t) = ai t, 
A(0) = Ai</.(t) 

Cl 



A20(O" 



(7.20) 



where (pi = ^(3 + 2oj), ai = y gq^i Ai = ^ and A2 = ci. This solution has an initial singularity, as we can 
see from the Ricci scalar 

3?= 12(2+0;) (7.21) 



X{(j)) increases with the time, in contradiction with the actual observations. However, for a particular combination 
of the today values of the free constants in this model we can get a small value of Aq 



-1 1/4 



to 



(7.22) 



J 



The today values of the corresponding deceleration and Hubble parameters are given by 

qo = 0, Ho = ^, (7.23) 

10 



such that to = I/Hq ~ 15.05 Gy. 



VIII. FINAL REMARKS 



We have considered a Brans-Dickc scalar-tensor theory, obtaining Friedmann- Robertson- Walker cosmological models 
with time dependent cosmological constant. The time dependence occurs in a natural way. Two ansatz were proposed: 
at^™ = a, with a constant, and X{t) = Ai0(t)"i -|- \2<j){t)'^^, in order to get exact inflationary solutions of the field 
equations, with a general state equation p = (7 — l)p. Our set of exact solutions depend on the values of 7, k, m and 

OJ. 

We classify the exact solutions of each case which we deal, according to the values of the free constants of our 
model. For vacuum with k = and m = 1/2, we get a non-relevant solution (according to our goal), with A = 0. 

For A; = 0, TO = 2/3, wc get a singular solution for which the scale factor decreases with the time as a{t) ^ t^^^"^, in 
accelerated way, but its predicted age is a negative value, then we conclude that this model has not physical meaning 
today. 

Furthermore for vacuum, we get for a flat case, an extended inflationary solution with initial singularity, and an 
exponential inflationary solution without singularity. For a not flat case we get a coasting singular solution. In all 
this models, the values of to and Cl^ (usually called Ha) are similar to the actual accepted values. 

We obtain exact solutions for a general equation of state p = (7 — l)/9. In the flat case {k = 0) we get an extended 
inflationary solution with initial singularity. The expansion of this model occurs in a accelerated way, independently 
of the equation of state. The values of to, Qm and Q0 depend on the value of the undetermined constants 7, and uj. 

The solution of the non-flat case cannot be expressed in terms of the cosmological time, but in terms of the parameter 
r. In such models the initial singularities can be avoided for some values of 7. 

As examples of our general solutions, we calculate the models for a dust and stiff matter fluid. For a flat dust model 
we get a slow extended inflationary solution with acceleration and time decaying cosmological constant. This solution 
has an initial singularity and its estimated age is to ^ 2/i?o, which is too big according to actual known values. A 
bigger growth rate is required in order to have smaller values of to and to be consistent with actual observations. 

On the other hand, for a non-flat dust model we get a slowly expanding solution a{t) ~ i^, with non-constant 
acceleration or deceleration depending on the relation between our free constants k/{3 + 2uj) > or fc/(3 -|- 2u)) < 
respectively. The models are singular or non-singular depending on the free constants a, k and w. 

As another application of our solutions with a general state equation, we consider a stiff matter fluid for which 
7 = 2. The solution of the corresponding flat case is a slowly expanding model a{t) ~ f^/^, with acceleration, without 
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singularities and with a time decaying cosmological "constant": A ^ t^/^. The vahdity of this model is restricted to 
some values of t. The actual age of the Universe predicted by this model, to ^ l/2Ho ^ 7.5 Gy, with Hq ~ 65 ± 5 km 
s^^Mpc^^; clearly this model is not applicable today. 

In the non-flat case for a stiff matter fluid, we get cosmological models for which the scale factor grows with the time 
from a minimum radius, and the cosmological "constant" decreases with the time, but as in the dust model, we have 
a set of free constants, whose values should determine the characteristics of this model. Thus, for k/{3 + 2u)) > we 
have a non-constant accelerated model without singularities, while for k/{3 + 2uj) < 0, our solution is a non-constant 
accelerated model with initial singularities and which validity is restricted to some values of t. In both cases the 
actual predicted age of the Universe, depend on the numerical value of the free constants. 

We solve separately the case of a radiation fluid (7 = 4/3) and a false vacuum fluid (7 = 0). For the radiation fluid 
we get a coasting model a{t) t, with a decaying cosmological "constant" and with initial singularity, independently 
of the curvature. The value of to ^ 15.05 Gy, obtained from this model is in fair agreement with actual observations. 

For a false vacuum fluid, we obtain a set of solutions which we classify depending on the range of values which 
our free constants may take. For the flat case, we get two family of solutions, one of them, corresponds to an 
exponential inflationary model with acceleration and without singularities. Another set of solutions is a kind of 
power law inflationary model a{t) ~ (c — ty, where e depend on the "free" constants m and to, restricted by physical 
requirements. For this last solution the cosmological constant is a binomial function of t, which decreases under 
specific conditions on the free constants. 

Additionally, we get a solution for the not flat case of a false vacuum fluid. In such a case we get a coasting 
model a{t) ~ t, which has initial singularity and with cosmological "constant" which is a binomial function of t: 
A ~ Xit~^ + \2t^- Such a cosmological "constant" would increases with the time in contradiction with the actual 
accepted value of A. For a particular combination of the today values of the free constants in this model, it is possible 
to obtain small values of Aq. The actual age predicted by this model is to ^ 15.05 Gy, which is again in fair agreement 
with the actual accepted value. 

Then, as we can see from the description of our exact solutions, some of them have not physical meaning today, some 
others are restricted to be valid during a specific period. Most of them are valid from an initial singularity until today, 
predicting an inflationary epoch, a cosmological "constant" which decreases with the time and the today observed 
acceleration, as well as an actual age of the universe which is in reasonable agreement with the actual observations. 
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